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ABSTRACT 



We present nine newly observed transits of TrES-3, taken as part of a tran- 
sit timing program using the RISE instrument on the Liverpool Telescope. A 
Markov- Chain Monte-Carlo analysis was used to determine the planet-star radius 
ratio and inclination of the system, which were found to be RpjR* = 0.1664to;ooig 
andi = 81.73^ respectively, consistent with previous results. The central tran- 
sit times and uncertainties were also calculated, using a residual-permutation 
algorithm as an independent check on the errors. A re-analysis of eight previ- 
ously published TrES-3 light curves was conducted to determine the transit times 
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and uncertainties using consistent techniques. Whilst the transit times were not 
found to be in agreement with a linear ephemeris, giving \ 2 — 35.07 for 15 degrees 
of freedom, we interpret this to be the result of systematics in the light curves 
rather than a real transit timing variation. This is because the light curves that 
show the largest deviation from a constant period either have relatively little out- 
of-transit coverage, or have clear systematics. A new ephemeris was calculated 
using the transit times, and was found to be T c (0) = 2454632.62610 ± 0.00006 
HJD and P = 1.3061864 ± 0.0000005 days. The transit times were then used to 
place upper mass limits as a function of the period ratio of a potential perturbing 
planet, showing that our data are sufficiently sensitive to have probed for sub- 
Earth mass planets in both interior and exterior 2:1 resonances, assuming the 
additional planet is in an initially circular orbit. 

Subject headings: methods: data analysis, stars: individual (TrES-3) , stars: 
planetary systems, techniques: photometric 



1. Introduction 



Transit surveys of extrasolar planets have vastly improved our understanding of plan- 
etary systems in recent years, with a rapid increase in the number of new discoveries^. 
Transiting systems are particularly important, because, when coupled with radial velocity 
measurements, they allow the measurement of the mass, radius and density of the planet. 
While the majority of these systems are Hot Jupiters, neither ground-based transit nor ra- 
dial velocity surveys have reached the precision required to search for Earth-sized planets. 
However, Earth-sized planets may be found via high precision ground-based observations 
through the detection of Transit Timing Variations (TTV). 

A transiting planet will maintain a constant period whilst orbiting its parent star (ex- 
cluding tidal effects and general relativity), unless acted on by a third body. Measuring 
the central transit times allows us to detect perturbations in the period thus revealing 
the presence of anoth e r body in the system, which is the princ i ple of the TTV metho d 
jMiralda-Escudel I2OO2I : iHolman fc Murravl I2OO5I : kgol et all I2OO5I : iHevl fc Gladmanl l2007h 
TTV is particularly sen sitive to small bodies in resonant orbits, or e ven exomoons (IKipping 
20091 ) and trojans (e.g. iFord fc Gaudil 120061 ; iFord fc Holmanl 120071 ). and therefore has the 



potential to provide the first detection of an Earth-sized body orbiting a main-sequence star 



http: / / exoplanet.eu/ 
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other than our own. 



Constraining a TTV signal requires many high precision light curves with high cadence. 
In theory we can measure TTVs to several seconds. However, in practice we are limited 
by correlated noise in the light curves, which may arise due to effects such as pixel-to- 
pixel sensitivity y ariations, temperature fluctuations, or changes in the observing conditions 
( jPont et al.l 120061 ). Indeed, there may also be non-instrumental effects caused by brightness 
variations of either the target or comparison stars. We are unable to distinguish any stellar 
activity from other sources of correlated noise in the light curves, and therefore they have 
the same detrimental effects as instrumental systematics when making transit observations. 

Providing that these sources of correlated noise can be kept to a minimum, it is still 
possible to measure central transit times to better than 10s (see §H}. This allows us to probe 
for the presence of Earth-sized planets in lo w-order mean-motion resonance, or more massive 
pertu rbers in non-resonant orbits (see e.g. ISteffen fc Agolll2005l ; lAgol fc Steffenl 120071 ; [Bean 
20091 ). 



RISE (Rapid Imager to Search for Exoplanets) is a fast camera mounted on the Liverpool 
Telescope (LT) on La Palma, primarily to obtain hig h precision light curves of transiting 
exoplanets (see ISteele et al.l 120081 ; iGibson et al.ll2008l ). It was commissioned in February 
2008, and observations of several exoplanet systems have been ongoing since then in an 
effort to detect TTV signals in these systems. 



Sozzetti et al.l (120091 ) presented eight transits of TrES-3, a G-type dwarf hosting a 1.9 M 



planet in a 1.3 day period (jO'Donovan et al.ll2007l l. They concluded that a linear period did 
not provide a particularly good fit to the central transit times, which indicates that either 
they underestimated the sytematics in their light curves, or that there is indeed a real TTV 
indicating a third body in the system. 

Here, we present a fur ther nine RISE transit light curves of TrES-3, and re- analyse those 
from lSozzetti et al.l (120091 ) using consistent techniques, in an effort to detect and understand 
any TTV signal. In §2] we describe the observations and data reduction, and in §3] describe 
how the light curves are modelled and in particular how the central transit times and uncer- 
tainties are found. Our results are presented in £JU and we use the transit timing residuals to 
place upper mass limits on a perturbing planet that could be present in the TrES-3 system 
without being detected from our observations. Finally, in §|5] we summarise and discuss our 
results. 
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2. Observations and data reduction 
2.1. RISE photometry 
Seven full and two partial transits were observed using the LT a nd RISE from 2008 



Mar ch 8 to 2008 August 4. The RISE instrument is described in detail in lSteele et al.l (120081 ) 



and iGibson et al.l (120081 ) . It consists of a frame transfer CCD which allows for continuous 
observation with effectively no dead time, a relatively large field of view (9.4 x 9.4 arcmin 
squared), and a single wide band filter (~500-700 nm). 

For all observations, an exposure time of 8s was used with the instrument in 2 x 2 
binning mode, giving a scale of 1.1 arcsec/pixel. For the full transits, 1350 images were 
obtained resulting in 3 hours of continuous observations, allowing ~50mins of observations 
both before and after the transit event. The images have a typical FWHM of ~2-4 pixels 
(~2.2-4.4 arcsec). The nights were clear for the majority of the observations, except for part 
of the nights of 2008 July 5 and 2008 August 4, where large scatter due to thin clouds can be 
seen towards the end of the light curves. A summary of the observations is given in Table [TJ 

Images were first debiased and flat fielded with combined twilight flats using standard 
IRAF0 routines. Aperture photometry was then performed on the target star and nearby 
companion stars using Pyrao and the DAOPHOT package. In each night different aperture 
sizes and numbers of comparison stars were used to minimise the out-of-transit RMS. These 
varied as the conditions and field orientation changed for each night of observations. 

The flux of TrES-3 was then divided by the sum of the flux from the companion stars (all 
checked to be non-variable) to obtain each lightcurve. Initial estimates of the photometric 
errors were calculated using the aperture electron flux, sky and read noise. The lightcurves 
were then normalised by dividing through with a linear function of time fitted to the out-of- 
transit data, setting the unocculted flux of TRES-3 equal to 1. The light curves, along with 
their best fit models and residuals (see §3.11) . are shown in Figures [1] and [2j 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

3 Pyraf is a product of the Space Telescope Science Institute, which is operated by AURA for NASA. 
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3. Light curve modelling and analysis 



3.1. Determination of system parameters 

In order to determine the system parameters f r om th e transit light curves, a param 



eterised model was constructed as in I Gibson et al.l (120081 ). This used Kepler's Laws and 
assumed a circular orbit to calculate the normalised separation (z) of the planet and star 
centres as a function of time from the stellar mass and radius (M* and R*), the planetary 
mass and radius (M p and R p ), the orbital period and inclination ( P and i), and finally a 



central transit time for each lightcurve (To >n ). The analytic models of iMandel fc Agoll (120021 ) 
were then used to calculate the stellar flux occulted by the planet from the normalised sepa- 
ration and the planet /star radius ratio (p) assuming the quadratic limb darkening function 



/ = l-a(l-/i)-6(l-/i) 2 , 
n 

where I is the intensity, /x is the cosine of the angle between the line-of-sight and the normal 
to the stellar surface, and a and b are the linear and quadratic limb darkening coefficients, 
respectively. 



Limb darkening parameters were obtained from the models of lClaretl (120001 ). We linearly 
interpolated th e ATLAS tables for T eff = 5 650K, logg = 4.4, [Fe/H] = -0.19 and v t = 
2.0 kms -1 (from lSozzetti et al.ll2009i ) to obtain limb darkening parameters in both the V and 
R bands. The average from the V and R bands was then adopted as our theoretical limb 
darkening parameters. Several tests were performed to examine the effects of the choice of 
limb darkening parameters on the results, which are described at the end of this section. 

A Markov-Chain Monte-Carlo (MCMC) alg orithm was then used to obtain the best 



fit parameters and their uncertainties (see e.g . , iTegmark et al.l 12004 ; iHolman et al.l 12006 



Winn et al.l 120081 ; ICollier Cameron et al.l I2QQ7I ) . This consists of calculating the \ 2 fitting 



statistic, 



JV 



X 



E 



{hobs - f h calc? , (M*-M ) 



+ 



where fj )0 bs is the flux observed at time j, Oj is the corresponding uncertainty and fj jCa ic is 
the flux calculated from the model for time j and for the set of physical parameters described 
above. The second term represents a Gaussi an prior placed on M +, where M and o"m are 
the stellar mass and uncertainty as given in ISozzetti et al.l (120091 ) . This allows the stellar 
mass to vary within constraints for each model fit, so that errors in the stellar mass are 
taken into account when extracting errors from the MCMC distributions. The stellar radius 

1/3 

was updated for each choice of M* using the scaling relation R+ oc MJ , whilst P and M p 



- 6- 



were held fixed at their previously determined values, as their uncertainties do not have any 
significant effect on the output probability distributions. Subsequent parameter sets are then 
chosen by perturbing small amounts to the previously accepted parameter set and are then 
accepted with probability exp(— A% 2 /2) at each point in the chain, where A% 2 represents 
the difference in \ 2 calculated for the o ld and new parameter sets. The procedure is the 
same as that used in I Gibson et al.l (120081 ). to which the reader is referred for details. 



To obtain reliable estimates of parameters and their uncertainties, it is important that 
the photometric errors are calculated accurately. The photometric errors <jj are first rescaled 
so that the best fitting model for each lightcurve has a redu ced y 2 of 1. It is also vital to 
accou nt for any correlated ("red") noise i n the data (see e.g., iPont et al.ll2006l ; iGillon et al. 



20061 ). The same procedure was used as in lGibson et al.l (120081 ). where we evaluat ed the pres- 
ence o f red noise in each light curve by calculating a factor (3 (> 1) according to lWinn et al. 



( 120081 ) and rescaled the photometric errors by this value. A value for (3 is determined by 
analysing the residuals from the best fit model of the lightcurves. Calculating the standard 
deviation of the residuals a±, and the standard deviation after binning the residuals into M 
bins of N points a at, one would expect 



0"! 



M 



y/N\ M~ 1 



in the absence of red noise. This is usually larger by a factor (3. However, the value deter- 
mined for (3 depends strongly on the choice of averaging time, i.e. M and N. Previously 
we have used an average of (3 values in the range 10-35 mins (the approximate time-scale 
of ingress or egress) to rescale the photometic errors. However, for this analysis we decided 
to use the maximum value for (3 in this range in order to be as conservative as possible in 
determining our resultant errors. 

Normalisation plays an important role in determining parameters and errors from light 
curves, and to account for this a further 2 parameters were added to the model for each 
transit. These were the out-of-transit flux (f 00 t,n) an d a time gradient (ta r ad,n), which are 
vital for TTV measurements as these affect the symmetry of the light curve and therefore the 
central transit times. An airmass correction was not used as previous studies have shown this 
produces similar res ults for full transits, but impedes chain convergence for partial transits 
^Gibson et allboosh . 



An initial MCMC analysis was used to estimate the starting parameters and jump 
functions for p, i, T 0n , f oot , n and tcrad,n- An MCMC run was then started for all nine 
light curves, fixing the central transit times to those determined in the initial run. Other 
parameters, such as the normalisation parameters, that are independent for each light curve 
were still allowed to vary. Five separate chains with 200 000 points were then computed 
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with the initial free parameters set by adding a 5cr gaussian random to their previously 
determined best fit values. The first 20% of each chain was eliminated to keep the initial 
conditions from influencing the results, and the remaining parts of the chains were merged 
to obtain the best fit values and uncertainties for each free parameter. The best fit value 
was set as the modal value of the probability distribution, and the la limits to the values 
where the integrals of the distribution from the minimum and maximum values were equal 
to 0.159. To test that the chai ns had all converged to the same region of parameter space, 
the Gelman & Rubin statistic (jGelman fc Rubin! Il992l ) was then calculated for each of the 
free parameters, and was found to be less than 0.5% from unity for all parameters, a good 
sign of mixing and convergence. 

To check for any errors that may have resulted from a poor choice of limb darkening 
parameters, the above procedure was repeated, this time allowing the linear limb darkening 
parameter (a) to vary freely whilst holding the quadratic limb darkening parameter (b) fixed 
at the theoretical value, as in lSouthworthl (120081 ) . This, ho wever, results in unph ysical models 
of the limb darkeni ng (negative values) , as was found in ISozzetti et al.l (120091 ) , but not by 
Gibson et al.l (120081 ) using the same technique for the WASP-3 system. This is probably due 
to the higher impact parameter of TRES-3, and therefore a higher sensitivity to the limb 
darkening parameters. As a com promise, the same a p riori constraint was imposed on the 
linear limb darkening as used in ISozzetti et al.l (120091 ). assuming that the limb darkenin g 
parameters do not drift from their theoretical values by more than 0.2 (jSouthworthl 120081 ) . 
This involved adding another term to the \ 2 function as follows; 



N 



X' 



I, Irf | (M,-M ) 2 | 



a , 



0.2 



where ao represents the theoretical limb darkening coefficient. As this causes significant 
increases in the errors determined for % and p, these results were adopted as our final system 
parameters. 

A further two checks were performed to test the limb darkening parameters. The first 
involved repeating this process by replacing the quadratic limb darkening coefficient deter- 
mined for the combined V+R filter by that obtained for the individual V and R filters. This 
causes no significant changes to our results. The second check involved allowing each light 
curve to have its own independently varying (linear) limb darkening coefficient (within the 
prior constraints), rather than having one set of limb darkening coefficients to describe all 
the light curves. This again caused no significant changes to our results. It is therefore 
favourable to have the same set of limb darkening coefficients to describe all of the light 
curves, as forcing the same transit shape may reveal systematics through small differences 
in each light curve, that could otherwise be hidden through varying the limb darkening 
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coefficients independently. 



3.2. Central Transit Times 

In order to calculate the central transit times for a TTV analysis we used the MCMC 
code, as described above, on each individual light curve, this time keeping the system pa- 
rameters p and i fixed at the best fit values determined in the previous section. Modeling 
the light curves individually has the advantage of needing much shorter chains, and does not 
result in underestimated uncertainties. This is because the central transit times are not very 
sensitive to the physical system parameters, but rather to those parameters that effect the 
symmetry of the light curves, in particular the norm alisation function. The same analysis 



was done on the light cuves from ISozzetti et al.l (120091 ). so that the central transit times and 



errors were found using consistent methods. 

For each light curve, five chains of length 50 000 were computed, and T and its uncer- 
tainty were extracted as before from the probability distribution after merging the chains 
(again discarding the first 20% of each). The linear limb darkening coefficient, stellar mass, 
and stellar radius were allowed to vary within the same constraints outlined before. Again, 
the Gelman & Rubin statistic was used to check for convergence. The systematics were 
accounted for using the same technique as described in the previous section, by re-scaling 
the errors of each light curve by a factor f3. Values for f3 are given in Tables [3] and HI 



A residual-pe rmutation (RP) or "prayer bead" algorithm (see e.g. ISouthworthl 12008 



Gillon et al.l 120091 ) was also used on each of the light curves to determine the errors in the 
transit times. This is another method commonly used to evaluate the effects of systematic 
noise on light curves, which often results in larger uncertainties than the MCMC method. 

The RP method consists of reconstructing the light curve by adding the residuals to the 
best fit model from the MCMC fit, each time shifting the residuals by a random amount, 
and performing a new fit on the light curve. 10 000 such fits were performed for each 
transit, with M*, R±, i and p selected from a Gau ssian distribut i on at the start of each 



using the stellar parameters and uncertainties from ISozzetti et al.l (120091 ). and the system 
parameters determined for the combined RISE light curves. The transit times, normalisation 
parameters, and linear limb darkening co-efficient were allowed to vary freely, using starting 
points determined randomly within 10<t from the best fit values. Errors in the central transit 
times were then estimated from the resulting distribution of fits. 



This method has the advantage that it preserves the actual correlated noise from the 
light curve, whereas the error re-scaling technique used alongside the MCMC fitting is sen- 
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sitive to the choice of averaging time. However, a comparison of the two methods showed, 
that in most cases, the errors from the MCMC fit were larger than those from RP. This 
is likely due to choosing the maximum value for f3 in the 10-35 min range to re-scale the 
photometric errors prior to the MCMC runs, rather than using the average value. 

For each transit the "worst case" was assumed, i.e. we adopted the error from the RP 
method only when it produced a larger uncertainty than the MCMC code. Note we always 
used best fit values from the MCMC fit as the RP method already assumed these transit 
times when reconstructing each light curve. The methods used to determine each of the 
timing errors are given in Tables [3] and HI 



4. Results 



4.1. System parameters 



The system parameters d erived from the MCMC fits of the RISE transits are shown in 
Table [2J ISozzetti et al.l (120091 ) undertook a thorough analysis of the stellar properties and 
radial velocities, and therefore we focus only on the planet parameters that are observable 
in the light curves, namely the inclination of the orbit and the ratio of the planet to stellar 
radius. We found i = 81.73lo!o4> an d P — 0.1664lo;ooi8- These are consistent with previously 
determined values, although with slightly smaller uncertainties, and theref ore the planet 
radiu s and density (derived from the stellar radius and planetary mass from ISozzetti et al. 
20091 ) are also consistent with previous studies. 



4.2. Transit ephemeris 

The central transit ti mes are shown i n Tab le [3] for the RISE light curves, and in Ta- 



ble H] for the light curves of ISozzetti et al.l (120091 ). where the transit times were found to be 



consistent to within ~ 0.3cr in all cases, and typically to less than O.lcr. Due to the more 
rigourous approach used to account for red noise in our analysis, the error bars were found 
to be ~ 10 - 40% larger. 

A new ephemeris was calculated by minimising \ 2 through fitting a linear function of 
Epoch E and Period P to the transit times 

T C (E) = T c (0) + EP, 

where E = was set to the transit from 2008 June 14 taken with RISE, as it has the smallest 
uncertainty. The results were T c (0) = 2454632.62610 ± 0.00006 and P = 1.3061864 ± 
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0.000 0005. Figure [3] shows a plot of the timing residuals of the RISE and ISozzetti et al. 
( 120091 ) transits using this updated ephemeris. 

For the RISE data a straight line fit yeilds x 2 = 13.49 for 7 degrees of freedom, and for 
the lSozzetti et al.l (120091 ) data gives x 2 — 19.40 for 6 degrees of freedom, much lower than the 
value of 35.22 found from their analysis, simply because of the larger timing errors. For the 
combined data set, x 2 — 35.07 for 15 degrees of freedom, and t herefo re a reduced x 2 of 2.34. 
This all seems to support the conclusions from ISozzetti et al.l (120091 ). that the uncertainties 
are underestimated, or that a linear period is not a good fit to the data. As we have been 
as sceptical as possible regarding the timing errors, this seems to suggest tentative evidence 
of a third body in the system perturbing the orbit of TrES-3b. 

However, after closer inspection of the light curves that contribute most to x 2 i this 
conclusion is less convincing. Most of the large contributors to x 2 are partial transits or 
have very little out-of-transit data. If we remove all of the light curves with less than 20 
mins of out-of-transit data either before ingress or after egress (transits E = -332, -319, -29, 
23 and 32) this results in a x 2 °f 13.53 for 10 degrees of freedom, or a reduced x 2 °f 1-35. 
This is because these transits are much more difficult to normalize due to lack of out-of- 
transit data, and unseen systematics could certainly cause the normalization gradient to be 
skewed, therefore effecting the symmetry of the light curves and hence the central transit 
times. This seems to suggest that transits need at least ~ 20 — 30 mins of out-of-transit data 
either side of the transit to be useful for transit timing studies, unless a more robust method 
of normalizing light curves and accounting for the errors is found. The largest remaining 
contributor to x 2 is transit E = 29 which lies ~ 2.4cr from the straight line fit. This transit 
not only has a high level of red-noise {(3 > 2), but a dip in the residuals from the best 
fit model is clearly seen around egress. The net effect of this on the model fit would be 
to "drag" the measurement of the central transit time later, as seen in the transit timing 
residuals. Removing this transit results in a reduced x 2 < 1> which suggests a constant 
period. Conclusions supporting a third body in this system would therefore rely on transits 
with little out-of-transit coverage and/or those with large visible systematics. 



4.3. Limits on a second planet in the TrES-3 system. 

Despite not revealing a significant TTV signal, the data can still be used to place upper 
mass limits on the presence of a hypothetical second planet in the TrES-3 system. The shape 
and amplitude of transit timing residuals are dependent on a large number of parameters, 
such as mass, period, eccentricity and argument of periastron of the perturbing planet. 
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In order to compute model timing residuals, the equations of motion for a three body 
system were integrated using a 4th order Runge-Kutta method. The first two bodies were 
set to represent the star and planet of the TrES-3 system, which was assumed to have an 
initially circular orbit. Transit times were then extracted when the star and the transiting 
planet were aligned along the direction of observation, with the third body representing the 
perturbing planet. The transit times were then fit with a linear function of time and the 
timing residuals used for comparison with the data. The orbits of the planets and direction 
of observation were assumed to be coplanar. 

Ideally we would like to search the parameter space of the perturbing planet completely 
and set upper mass limits at each point. However, this is not possible given the large amount 
of computation required to produce a model of timing residuals at each point in such a large 
parameter space. Therefore, some simplifications and assumptions were made. Firstly, we as- 
sumed that the amplitude of the timing residuals for a g iven perturbing orbital configuration 



are p roportional to the mass of the perturbing planet (lAgol et al.ll2005l ; iHolman fc Murray 



20051 ). and verified this by constructing models with a range of perturbing planet masses. 
Secondly, we assumed that the perturbing planet had a starting eccentricity of 0, as an in- 
crease in eccentricity generally increases the amplitude of the timing residuals and therefore 
to set upper masses as a function of period we only need to investigate perturbing planets 
on circular orbits. This assumption is tested later in this section. 

Models were created for an Earth-massed perturbing planet with a period ratio dis- 
tributed from 0.2 to 5.0 (the regime in which relatively small masses may be detected), 
increasing the sampling around both the interior and exterior 2:1 resonances, where we ex- 
pect to probe for the smallest masses. For each model produced, the transit times were 
extracted for a range of observation directions. 

To calculate the maximum allowed mass for each model, x 2 was calculated by fitting 
the model residuals to the measured timing residuals from the light curves. The mass of 
the perturbing planet was increase d (or decreased) by s c aling the timing resi duals until x 2 



was increased by a value A% = 9 (ISteffen fc Agoll 120051 ; lAgol &: Steffenl 120071 ) from that of 



a constant period (ie timing residuals = 0) which corresponds to a 3<r confidence limit. We 
then minimised x 2 along epoch only, and then let the mass of the perturber grow again until 
the maximum allowed mass for each model was determined. This procedure was repeated 
for the range of observation directions of each model, and the largest upper mass determined 
was assumed as our upper mass limit for each period ratio. 

To check that the starting mass of each model had no impact on the mass limits found 
for each period (i.e. test that residuals are indeed proportional to the perturbing mass), 
models were re-calculated with the mass of the perturbing planet set as the upper mass 
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limits found from the \ 2 fits. Upper masses limits were then determined as before. This 
process was repeated twice and was found to make little difference to the final upper mass 
limits, therefore justifying our assumption. 

Figure H] shows a plot of the resulting upper mass limits found as a function of the period 
ratio. The solid black line shows the upper mass limits found for the 3 body simulations, and 
the horizontal dashed line represents an Earth-mass planet. The results show that we have 
probed for masses as low as 0.97 M e and 0.71 M e in the interior and exterior 2:1 resonances, 
respectively. 

To test our assumption that perturbers on initially circular orbits will cause the small- 
est perturbations, and therefore can be used to set upper mass limits on a perturbing mass 
as a function of period ratio, we created models this time allowing the perturbing mass to 
have non-zero initial eccentricity. A set of models was created with a period range spanning 
the exterior 2:1 resonance and the perturbing bodies eccentricity ranging from to 0.15. 
It was found that generally the amplitude of the signal drops and reaches a minimum be- 
tween e ~0.01 and 0.12, before increasing again, and could drop by as much as an order 
of magnitude. A similar (but smaller) effect was found for the interior 2:1 resonance. This 
invalidates our assumption, and suggests that to set upper mass limits around resonance at 
least the period, eccentricity and argument of periastron of the perturbing planet needs to 
be explored in parameter space, which is beyond the scope of this paper. However, upper 
mass limits were estimated using this set of models and the same technique as before, and 
we found that more realistic upper mass limits are ~ 3 — 4M e and ~ 10 — 15M e in the 
interior and exterior 2:1 resonances, respectively. Out of resonance, it was found that the 
amplitude of TTV signals increases with eccentricity of the perturbing planet, and thus our 
assumption and upper mass limits are valid. 



A long term stability analysis was not performed for the 3-body systems. iBeanl (120091 ) 
found that for the CoRoT-1 system, only test particles with period ratios greater than ~ 1.8 
were stable for more than 10 6 orbits of the transiting planet. CoRoT-1 has a similar G-dwarf 
host star and a slightly longer period (~ 1.5 days) , which suggests that a similar analysis 



would prove useful here. iBarnes fc Greenbergi (120061 ) explore the stability limits in exoplanet 



systems, and provide an inequality to test whether a system is Hill stable (equation 2). Using 
this inequality for the TrES-3 system (assuming an Earth-massed perturber), places lower 
and upper limits on the period ratio of 0.64 and 1.59, respectively. The resulting region 
not guaranteed to be Hill stable is marked on figure H] by the grey shading, although stable 
configurations may still occur in this region. Trojan compan ions could also exist in stable 



orbits near the 1:1 resonance. iMadhusudhan &: Winn! (120091 ) placed a 2a upper mass limit 



of 81.3 Earth masses on a trojan in the TrES-3 system by combining transit observations 
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and radial velocity data. 

5. Summary and discussion 

This paper presents the first transits taken using RISE specifically for a transit timing 
analysis, consisting of nine light curves of TrES-3. The transits were fit with an MCMC 
code and the derived system parameters found to be consistent with previous studies. Two 
different methods were used to determine the errors in the central transit times, trying to 
take into account the systematics in the light curves. These were scaling the errors bars 
by a constant prior to MCMC fitting after analysing the residuals, and using a residual 
permutation algorithm. The largest error found was used for each transit. We have shown 
that when systematics are kept at a minimum, it is possible to determine transit times to 
~10 seconds - the level of accuracy expected from RISE. 

Whilst the transit times appear to deviate significantly from a constant period when 
a x 2 analysis is performed, those that contribute most to the deviations tend to have very 
little out of transit data or obvious systematics. After removing transits with less than ~20 
mins of out-of-transit coverage either before ingress or after egress, the data are consistent 
with a constant period, and therefore no conclusive evidence was found for the presence of a 
second planet in the TrES-3 system. The transit times were then used to place upper mass 
limits on a perturbing planet as a function of period that could be present in the system yet 
not detected through our observations. This showed that our observations were sensitive to 
Earth- mass planets or smaller in the interior and exterior 2:1 resonances when we assume 
the perturbing mass is on an initially circular orbit. However, larger planets may exist in 
low eccentricity orbits around the 2:1 resonances, and exploring period, eccentricity and 
argument of periastron in parameter space is called for to set true upper mass limits, as well 
as a long term stability analysis. 

This study highlights the difficulties in attempting to detect a TTV signal, as we need 
to have complete confidence in the error bars calculated, which may require a more robust 
method to normalize the light curves and deal with red noise, especially if we are to trust 
transits with limited out-of-transit data. Confirming a true TTV signal may therefore re- 
quire some obvious structure in the observed residuals (which we would expect for resonant 
systems), rather than relying on a larger than expected "scatter" of points. 

RISE was designed and built with resources made available from Queens University 
Belfast, Liverpool John Moores University and the University of Manchester. The Liverpool 
Telescope is operated on the island of La Palma by Liverpool John Moores University in the 
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Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias 
with financial support from the UK Science and Technology Facilities Council. D.L.P. was 
supported by a Leverhulme Research Fellowship for the duration of this work. F.P.K. is 
grateful to AWE Aldermaston for the award of a William Penney Fellowship. We also thank 
A. Sozzetti, for making his data available for re-analysis, and the referee, D. Fabrycky, for 
comments which improved the content of this paper. 

Facilities: Liverpool: 2m (RISE) 
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Fig. 1.— RISE light curves of TrES-3 taken from 2008 March 8 to 2008 July 5 with their 
best fit models from the MCMC analysis over-plotted. Residuals from the best fit model are 
shown below each light curve offset. 
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Fig. 2. — Same as Figured], for light curves from 2008 July 14 to 2008 August 4. 
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Fig. 3. — Timing residuals of the RISE transits (triangles) and those from ISozzetti et al 
(120091 . squares). 
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Fig. 4. — Upper mass limits of a hypothetical 2nd planet in the TrES-3 system as a function of 
period ratio. The solid black line represents the upper mass found for the 3 body simulations, 
and the horizontal dashed line represents an Earth-mass planet. The region where an Earth- 
massed planet is not guaranteed to be Hill stable is marked by the grey shading. 
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Table 1: Summary of the RISE light curves of TrES-3. 



Night 


No. exposures 


No. comparison 


Aperture size 


RMS (residuals) 






stars 


(pixels) 


(mmag) 


2008 Mar 8 


1350 


4 


8 


1.31 


2008 May 28 


1350 


4 


7 


1.42 


2008 Jun 14 


1350 


7 


5 


1.00 


2008 Jul 1 


1350 


6 


4 


1.10 


2008 Jul 5 


1350 


7 


6 


3.31 


2008 Jul 14 


825 


6 


7 


1.81 


2008 Jul 22 


1350 


2 


4 


1.57 


2008 Jul 26 


1125 


8 


4 


1.18 


2008 Aug 4 


1350 


2 


4 


2.60 



Table 2: Parameters and la uncertainties for TrES-3 as derived from MCMC fitting of RISE 
light curves and some further calculated parameters. 



Parameter 


Symbol 


Value 


Units 


Planet /Star radius ratio 


P 


0-16641°;°^ 




Orbital inclination 


i 


81.73ig^ 


deg 


Impact parameter 


b 


n oc:o+0.004 

u.odz_ 013 




Transit duration 


T d 


1 QQO+0.024 

i - ooz -o.oio 


hours 


Transit epoch 


To 


2454632.62610 ± 0.00006 


HJD 


Period 


P 


1.3061864 ±0.0000005 


days 


Planet radius 


R p 


•. OA 1+0.025 
—0.035 


Rj 


Planet mass a 


M p 


l qin +a075 

- L - y - LU -0.080 


Mj 


Planet density 


P P 


n voo+o- 047 

U. 1 J^_ .042 


Pj 


Planetary surface gravity 


log 9 P 


O /lOI +0.023 
O.^X_ 002 2 


[cgs] 



a From Sozzctt i et al. 



(|2009l ). displayed here for convenience. 
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Table 3: Central transit times and uncertainties for the RISE photometry including the error 
source. 



Epoch 


Central Transit Time 


Uncertaintly 


P h 


Error Source 




[HJD] 


(days) 






-75 


2454534.66243 


0.00017 


1.07 


RP 


-13 


2454615.64553 


0.00017 


1.49 


MCMC 





2454632.62613 


0.00011 


1.25 


MCMC 


13 


2454649.60634 


0.00013 


1.31 


MCMC 


16 


2454653.52504 


0.00037 


1.95 


MCMC 


23 


2454662.66896 


0.00034 


1.56 


MCMC 


29 


2454670.50630 


0.00033 


2.04 


RP 


32 


2454674.42423 


0.00052 


2.83 


RP 


39 


2454683.56734 


0.00018 


1.00 


MCMC 



^Re-scale factor from red noise analysis (See $3J}- 



Table 4: Central transit times and uncertainties for the light curves of ISozzetti et al.l (120091 ) 
including the error source. 



Epoch 


Central Transit Time 


Uncertaintly 


(3 C 


Error Source 




[HJD] 


(days) 






-342 


2454185.91040 


0.00028 


1.63 


MCMC 


-332 


2454198.97307 


0.00033 


1.52 


MCMC 


-320 


2454214.64631 


0.00036 


1.58 


MCMC 


-319 


2454215.95210 


0.00024 


1.29 


MCMC 


-74 


2454535.96825 


0.00023 


1.39 


MCMC 


-61 


2454552.94898 


0.00021 


1.52 


MCMC 


-48 


2454569.92910 


0.00021 


1.34 


MCMC 


-29 


2454594.74597 


0.00028 


1.30 


MCMC 



c See Table El 



